Dynamics of Pinned Magnetic Vortices 
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We observe the dynamics of a single magnetic vortex in the presence of a random array of pinning 
sites. At low excitation amplitudes, the vortex core gyrates about its equilibrium position with a 
frequency that is characteristic of a single pinning site. At high amplitudes, the frequency of gyration 
is determined by the magnetostatic energy of the entire vortex, which is confined in a micron-scale 
disk. We observe a sharp transition between these two amplitude regimes that is due to depinning 
of the vortex core from a local defect. The distribution of pinning sites is determined by mapping 
fluctuations in the frequency as the vortex core is displaced by a static in-plane magnetic field. 

PACS numbers: 75.75.+a,75.40.Gb 
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The excitation spectrum of a sub-micron ferromagnetic 
particle is influenced profoundly by its shape. An exam- 
ple of considerable interest is the vortex state of a soft 
ferromagnetic disk, in which the lowest frequency exci- 
tation is a translational mode of the vortex core. This 
mode, in which the vortex core gyrates about its equilib- 
rium position, has been studied in several recent experi- 
ments P, 0, El 0, H, El ■ Although the core has a diameter 
of the order of the exchange length (~10 nm) 0,0], the 
eigenfrequency of the gyrotropic mode is determined only 
by the geometrical aspect ratio (diameter over thickness) 
of the disk in which it is confined In particular, 

the gyrotropic frequency should be independent of the 
location of the vortex core in the disk |fj. As demon- 
strated recently, however, the vortex core can be pinned 
by defects in a patterned thin film 0, 0, 0] . In the 
case of intrinsic defects, the range of the effective pinning 
potential appears to be on the order of ten nanometers 
hj. This raises the question of how a nanoscale defect 
influences the vortex gyrotropic motion. 

In this Letter, we report that pinning has a pronounced 
effect on magnetic vortex dynamics at small amplitudes. 
The vortex gyrotropic frequency in 50 nm thick permal- 
loy (Py) disks with diameters from 600 nm to 2 //m fluc- 
tuates by a factor of at least two as the core is displaced 
by a static magnetic field over length scales ~ 10 nm. 
The dependence on core location indicates that the small- 
amplitude dynamics are influenced strongly by the char- 
acteristics of a particular pinning site. Each site has a 
critical excitation amplitude above which the vortex be- 
comes depinned and gyrates at the frequency determined 
by the aspect ratio of the disk. We image the spatial 
distribution of pinning sites in a disk by displacing the 
vortex core in two dimensions while monitoring the gy- 
rotropic frequency. Using this approach, we find a density 
of up to ~ 2 x 10 11 cm~ 2 for the pinning sites in our films. 
Although this density cannot be correlated directly with 
the physical morphology of the permalloy films, the ex- 
treme sensitivity to core displacement indicates that the 
magnetic microstructure fluctuates strongly on nanome- 



ter scales. 

The samples were patterned from 50 nm thick films 
of Py (Nio.8iFeo.19), which has negligible magnetocrys- 
talline anisotropy. To enhance pinning, crystallographi- 
cally textured films were sputtered at a base pressure of 
5 x 10 -8 Torr onto a buffer/ substrate system of 30 nm 
Cu(lll)/ Al 2 O3(0001), in which the grain size depends 
on substrate temperature ^4|. The films were capped 
at room temperature with 2.5 nm Al and 50 nm (op- 
tically transparent) SiN. Results will be presented for 
two samples. A small-grain (SG) sample was grown at 
30° C, yielding an average grain diameter of « 35 nm as 
determined by atomic force microscopy. A large-grain 
(LG) sample with grain diameter rj 85 nm was grown 
at 200°C. The average roughness (Ra)is « 4 A (SG) and 
w 8 A (LG). Wide angle X-ray diffraction indicates (111) 
texture and low strain. Patterning was performed by 
electron beam lithography and dry etching. 

We use time-resolved Kerr microscopy (TRKM) [lU to 
study vortex core dynamics in individual magnetic disks. 
The substrate is polished to a thickness less than 30 /im, 
and the sample is then positioned on a coplanar waveg- 
uide as shown in Fig. ^a). The sample is excited by a 
magnetic field pulse with a temporal width less than 120 
ps and amplitude 2 to 13 Oe (estimated from the mag- 
nitude of the current) oriented in the plane of the disk. 
A linearly polarized laser pulse (duration 150 fs, A = 810 
nm) is focused through an oil-immersion objective to a 
spot with a full width at half maximum (FWHM) ~ 400 
nm. Using the polar Kerr effect, we measure the pump- 
induced change in the z-component of the magnetization 
as a function of the time delay between the pump and 
the probe. In the vortex state, shown schematically in 
Fig. IHb), the torque due to the pump pulse rotates the 
magnetization in the upper and lower halves of the disk 
into and out of the plane. Fig. QJc) shows a polar Kerr 
image of a 2 /im diameter SG disk at the peak of the 
pump pulse, demonstrating the expected contrast. The 
left half of this image is data and the right half is a mi- 
cromagnetic simulation of the response convolved with 
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FIG. 1: (color online) (a) Schematic of the experiment, (b) 
Schematic of a vortex. The z-component (out-of-plane) of the 
magnetization is non-zero only in the core, (c) In-plane mag- 
netization configuration for a 2 /im disk (arrows) and image 
of the ^-component of the magnetization (shading) 80 psec af- 
ter the arrival of the field pulse. The split frame shows both 
experimental (left) and simulated (right)images in a field of 
Oe. (d) The same image in a static field of 100 Oe. (e) 
Time scans of the polar Kerr signal obtained in several dif- 
ferent magnetic fields with the objective positioned 300 nm 
from the center of a 1 /im diameter disk. 



a Gaussian of FWHM 400 nm [16j. An image taken in 
a field of 100 Oe applied in the ^-direction is shown in 
Fig. Old) . As expected, the vortex core is displaced in 
a direction perpendicular to the field 0, El ■ Measure- 
ments of the position y of the nodal line in the image 
allow for a calibration of the core displacement, which 
is approximately linear with field. Experimental calibra- 
tions are dy/dH — 4.5 ± 0.3 nm/Oe for the 2 /zm disk 
and 0.9 ± 0.2 nm/Oe for the 1 /im disk. In comparison, 
micromagnetic simulations give dy/dH — 3.9 nm/Oe and 
1.1 nm/Oe respectively. dy/dH increases approximately 
quadratically with disk diameter in both experiment and 
simulations, although we cannot measure the displace- 
ment reliably for the smallest disks. 

Time scans obtained at a radius of 300 nm from the 
center of a 1 /im diameter SG disk are shown in Fig.^e). 
The oscillations are due to the gyrotropic motion of the 
vortex core about its equilibrium position 0, 0, El 01 • 
The sub-GHz gyrotropic frequency ujq and damping times 
greater than 10 ns observed in Fig. E^e) are typical 
III L3- LI , but the dependence of the frequency on applied 
field is unexpected. In the case of an ideal vortex [9j, 
uj = 7M s £ 2 /2x, where 7 is the gyromagnetic ratio, M s 
is the saturation magnetization, £ is a constant reflecting 
the boundary conditions, and x = dM/dH is the sus- 
ceptibility, which is nearly independent of field 0] . The 
gyrotropic frequency should therefore show very weak de- 
pendence on magnetic field, as confirmed in a recent res- 
onance experiment 

We have also investigated the field dependence of ojq 
for different excitation amplitudes. The Fourier trans- 
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FIG. 2: (color online) Field dependence of the vortex gy- 
rotropic frequency for 1/im (a) and 2/im (b) diameter disks 
fabricated from SG films. The amplitudes of the excitation 
pulse are 5 Oe (black squares), 11 Oe (red circles), and 13 Oe 
(blue triangles) . The dotted lines are the expected frequencies 
from micromagnetic simulations. Inset: The average period 
of the fluctuations (black squares) vs disk diameter. The solid 
line is a model (see text) based on a defect spacing of 25 nm. 



form of the time-resolved data at each field are fit to 
a Lorentzian to obtain u)q. Fig. [2] shows the depen- 
dence of ujq on static field for 1/zm and 2/im diameter SG 
disks. Each set of data corresponds to a different pump 
pulse amplitude h y . For small pump pulse amplitudes, 
the observed gyrotropic frequency fluctuates irregularly, 
with frequency peaks that are twice the expected field- 
independent value (dotted lines in Fig. [2J 9, 10]. The 
fluctuations are highly repeatable for a given disk, but 
are uncorrelated between neighboring disks. The data 
are independent of the vortex polarity (orientation of the 
vortex core into or out of the plane) but the field depen- 
dence is inverted upon reversal of the chirality (curling 
direction of the in-plane magnetization). The frequency 
fluctuates with a characteristic period AH — 23 Oe (ob- 
tained by counting frequency peaks over a field range of 
hundreds of Oe) for the 1 /zm disk and 7 Oe for the 2 /im 
disk. AH measured for different disk diameters is shown 
as symbols in the inset of Fig. [5] Note that AH increases 
by approximately a factor of 7 as the disk diameter de- 
creases from 2 /im to 600 nm. 

As shown in Fig. [2 an excitation pulse amplitude of 
13 Oe slightly suppresses the field-dependent frequency 
fluctuations for the 1 /im disk and almost completely 
eliminates the fluctuations for the 2 /zm disk. The de- 
pendence of the magnitude of the response on pulse am- 
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FIG. 3: (color online) The vortex gyrotropic frequency (black 
squares, left axis) for the 2 /im disk is shown as a function 
of field pulse amplitude at a fixed static field of -23 Oe. The 
maximum magnitude of the gyrotropic response (blue circles, 
right axis) increases rapidly for pulse amplitudes above 8 Oe. 
Time-resolved data at several different excitation pulse am- 
plitudes are shown in the inset. 



plitude is shown for the 2 disk in Fig. [3J For this 
measurement the static field is held fixed at one of the 
frequency peaks of Fig. [S^b) while the excitation pulse 
amplitude is varied. The most important feature of Fig-El 
is the rapid drop in the gyrotropic frequency from 0.4 to 
0.2 GHz at a pulse amplitude of 8 Oe. The magnitude 
of the response increases rapidly above this critical pulse 
amplitude. 

The important features of Figs. |2l and [3] are consistent 
with the existence of point-like defects in the Py disk 
that pin the vortex core. As shown in a recent Lorentz 
microscopy experiment jl3l| . the core moves from defect 
to defect as the magnetic field is varied. Fig. indicates 
that movement of the core among pinning sites leads to 
fluctuations in the gyrotropic frequency. Different disks 
fabricated from the same film should have an identi- 
cal areal defect density n, and therefore the fluctuations 
should be associated with a single length scale. To test 
this hypothesis, we use the calibrations dy/dH for the 
core displacement determined from micromagnetic sim- 
ulations for each disk diameter. As noted above, these 
agree well with the experimental calibrations measured 
for the larger disks. For each diameter d, we determine 
the field scale AH = Ay(dy / 'dH)^ 1 required to trans- 
late the vortex core over a distance Ay = y/ (1/n). The 
results for Ay = 25 nm are shown as the solid curve in 
the inset of Fig. [5] The good agreement with experiment 
demonstrates that the different periods AH for the dif- 
ferent diameters are associated with a single length scale 
characteristic of the defect density in the Py film. 

The second important feature of the experimental data 
is the existence of a critical pulse amplitude above which 
the gyrotropic frequency decreases abruptly to the value 
expected given the magnetostatic potential of the disk 
0. The critical amplitude corresponds to the depinning 



of the vortex core from a defect, which occurs when the 
radius of the gyrotropic orbit exceeds the range of the 
pinning center. At the maximum pulse amplitude, we 
are able to achieve depinning for almost all sites in disks 
with diameters greater than 1.5 /im. Above the depin- 
ning transition, the dynamics should be governed by the 
Thiele force equation [HHOH 



in 



where G = — Gz is the vortex gyrovector, X is the dis- 
placement of the vortex core from its equilibrium posi- 
tion, k is a stiffness determined from the magnetostatic 
energy of the displaced vortex, and D is a damping con- 
stant. F(t) _L h(i) is the time dependent excitation force, 
with h(t) the excitation field. The resonant gyrotropic 
frequency obtained from the solution of Eq.^is u>o = k/G 
|l9|. and the maximum response amplitude is propor- 
tional to 1/wo-D- It has been shown elsewhere that the 
solutions of the Thiele equation provide an excellent de- 
scription of the gyrotropic mode in permalloy disks 0,0- 
For the 2 /j,m disk, we find lu is independent of static 
magnetic field for pulse amplitudes >10 Oe, which corre- 
sponds to a gyrotropic orbital radius > 10 nm. The ob- 
served field independence in this high amplitude regime 
is consistent with the results of Novosad el al.]jfy. 

Except for the higher frequency, the vortex gyrotropic 
motion in the pinned case appears to be qualitatively 
similar to the unpinned regime. We therefore consider 
the possibility that Eq.^also applies for pinned vortices, 
but with a higher effective spring constant k 20]. For 
the case shown in Fig. |3 the change in ui at the crit- 
ical pulse amplitude implies that k in the pinned case 
is larger by a factor of 2. Changes in k alone, however, 
are not sufficient to account for the rapid increase in the 
magnitude of the response above the critical pulse am- 
plitude. Furthermore, we do not observe a clear trend 
in the damping times for pinned and unpinned vortices, 
and it is unlikely that G, which is independent of the size 
and shape of the vortex core, can change sufficiently to 
account for the observed increase in the gyrotropic signal. 
We therefore do not have a complete explanation of the 
dramatic increase in the signal above the critical pulse 
amplitude. Interestingly, the signals at high amplitudes 
are noticeably non-sinusoidal. This is probably due to 
deformations in the orbit caused by nearby pinning sites 
and not by any intrinsic non-linearity. 

By applying orthogonal in- plane magnetic fields, it is 
possible to move the vortex core anywhere in the plane 
of the disk. The core thereby functions as a nanoscale 
scanning probe, and the local value of the gyrotropic fre- 
quency can be used to map the pinning potential. We 
have used this approach to scan 1 fim diameter SG and 
LG samples using a pulse amplitude of « 7 Oe. Typi- 
cal atomic force micrographs (AFM's) of the SG and LG 
films are shown in Figs.QIa) and (b). Magnetic fields in 
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FIG. 4: (color online) (a),(b) Representative atomic force mi- 
croscopy images of the SG (a) and LG (b) samples. (c),(d) 
Contour maps of the gyrotropic frequency as a function of 
static field applied along both in-plane directions for 1 /im 
diameter SG (c) and LG (d) samples respectively, obtained 
at a pulse amplitude of ~ 7 Oe. (e) Contour map of the gy- 
rotropic frequency for the SG sample with the vortex chirality 
reversed. The image is almost identical to (c) after rotation 
by 180°. (f) Scale bar for frequency in (c) - (e) along with a 
representation of the vortex core in a random array of defects. 



the x and y directions are varied in steps of 5 Oe over a 
range of 125 Oe, which leads to a core displacement over 
a spatial region 140 nm x 140 nm (the size of the AFM 
images). Maps of the gyrotropic frequency for 1 ^m disks 
of each film are shown in Figs.^c) and (d). The pinning 
sites appear as regions of high frequency. Changing the 
vortex chirality, as for the case of Fig.^Je), reverses the 
direction of core displacement for a given field. The im- 
age of Fig. 0Jc) is therefore identical to that of Fig. Hfb) 
after inversion about the origin. 

From the images of Fig.^J we infer the density of pin- 
ning defects to be m 1.7 x 10 11 cm~ 2 for the SG sample 
and ~ 1.0 x 10 11 cm~ 2 for the LG sample. The areal den- 
sity of grains in Figs. 0Ja) and (b) is rj 8 x 10 10 cm~ 2 (SG) 
and k 1.4 x 10 10 cm~ 2 (LG). Although the density of pin- 
ning sites decreases with increasing grain size, there does 
not appear to be a direct correlation between the grain 
structure and the location of pinning sites. One might 
expect the strongest pinning sites to be at grain bound- 



aries, but the regions of high frequency in Figs.^Jc) and 
(d) appear point-like rather than filamentary. Using the 
core displacement calibration dy/dH, the typical range 
of a pinning potential in the SG sample is approximately 
10 nm, which is of the order of the core diameter. The 
point-like contrast in Fig. 4 as well as the threshold be- 
havior of Fig.|21suggest that the pinning potential acts on 
the core itself, lowering its exchange and magnetostatic 
energy. We believe that the energy associated with the 
in-plane spins in the vicinity of the core is nearly inde- 
pendent of the core position, although small fluctuations 
of the in-plane anisotropy energy from grain to grain have 
also been proposed as a source of vortex pinning 0] . 

In summary, we have demonstrated that pinning asso- 
ciated with intrinsic defects has a profound influence on 
the dynamics of magnetic vortices. For small amplitudes, 
the field-dependence of the gyrotropic mode frequency is 
governed by the density and strength of pinning sites. In 
this manner the vortex core serves as a sensitive probe of 
magnetic microstructure. 
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